Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Phenomenological Theory of Multiple Spin Density Waves 
in fee Transition Metals 



O 
O 

(N 



(N 



-a 
c 

o 
o 



> 
o 

o 

O 

-i— > 



i 

c 

o 
o 



X 



Takashi Uchida 1 * and Yoshiro Kakehashi 2 

1 Hokkaido Institute of Technology, Maeda 7-15-4-1, Teine-ku, Sapporo 006-8585 
2 Department of Physics and Earth Sciences, Faculty of Science, University of Ryukyus, 1 Senbaru, Nishihara, Okinawa, 

903-0213 

(Received February 6, 2008) 

The relative stability among the multiple spin density wave (MSDW) states in fee transition 
metals has been investigated on the basis of a Ginzburg-Landau type of free energy with 
terms up to the fourth order in magnetic moments. Obtained magnetic phase diagrams in 
the space of expansion coefficients indicate the possibility of various 3Q MSDW states in fee 
transition metals: the commensurate 3Q state with Q = 2n/a (a being the lattice constant), 
the incommensurate linearly polarized 3Q state, and the incommensurate helically polarized 
3Q state. It is shown that these 3Q states are always stabilized, when they are compared with 
the corresponding 2Q and 1Q states, and their magnetic moment amplitudes are the largest 
among those of the three states. The results are compared with previous results of the ground- 
state calculations, and possible scenarios to explain the experimental data of cubic 7-Fe are 
proposed. 
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1. Introduction 

Transition metals with nearly half-filled bands, such as 
Cr, 1 Mn, 2 ~ 4 7-Fe, 5 ' 6 and their alloys, 7-12 show complex 
magnetic structures due to competing magnetic interac- 
tions. The determination of their magnetic structures has 
long been a challenging problem in both theory and ex- 
periment in the study of metallic magnetism. 13 Of these, 
iron in the fee phase (7-Fe) has received special atten- 
tion, since it is located at a crossover point from the 
ferromagnetic to the antiferromagnetic state on the pe- 
riodic table and has been suggested to show spin density 
wave (SDW) states. 

Early experimental data on 7-Fe precipitates on a 
Cu matrix 14 and the extrapolation from those on 7- 
FcMn alloys 9 suggested that bulk 7-Fe shows the first- 
kind antiferromagnetic (AF) structure with wave vector 
Q = (0, 0, l)2w/a. Here a denotes the lattice constant. 
Later, the neutron diffraction measurements 5 on cubic 7- 
Feioo-xCo x (a; < 4) alloy precipitates in Cu showed mag- 
netic satellite peaks at wave vector Q = (0.1,0, l)2%/a. 
The magnetic structure was suggested to be a helical 
SDW, but was not determined precisely because of the 
high symmetry of the crystal structure. On the other 
hand, thin Fe films epitaxially grown on Cu were re- 
ported to show simple ferromagnetism 15 ' 16 or the co- 
existence of a low-spin AF state and a high-spin ferro- 
magnetic state, 17-19 depending on the film thickness. Re- 
cent studies 6 on Fe films suggested, however, that the fee 
phase of Fe is formed only for film thicknesses of 5 to 11 
monolayers and SDW is realized in these films. 

Theoretically, the magnetic structure of bulk cubic 7- 
Fe was investigated intensively by means of the ground- 
state theories of electronic-structure calculations. In 
most of the calculations, 20-26 the IQ helical SDW struc- 



ture was assumed for the bulk cubic 7-Fe and the wave 
vector that minimizes the total energy was determined. 
The obtained wave vectors, however, are different among 
different theoretical approaches. The linear muffin-tin 
orbital (LMTO) calculations by Mryasov et al. 20 and 
the augmented spherical wave (ASW) calculations by 
Uhl et al., 21 both based on the local density approxi- 
mation (LDA), yielded ground-state wave vector Q = 
(0, 0, 0.6)27r/a for lattice constant a = 6.8 a.u. On the 
other hand, Korling and Ergon 22 performed the LMTO 
calculations using the generalized gradient approxima- 
tion (GGA), and found the energy minimum at Q = 
(0.5, 0, l)2%/a. Furthermore, the recent full potential cal- 
culations suggest the possibility of other ground-state 
wave vectors. Bylander and Kleinman 23-25 performed full 
potential calculations using the ultrasoft pseudopotential 
and found the energy minimum at Q = (0, 0, 0.55)27r/a. 
By means of a modified ASW method, Knopfle et al. 26 
found the energy minimum at Q — (0.15, 0, l)27r/a for 
lattice constant a < 6.75 a.u. Although the wave vector 
obtained in their calculations is close to the experimental 
value, more recent ground-state calculations 27 based on 
the full potential augmented-plane-wave method, which 
accounts for more complex magnetic structures, show 
that the 1Q helical SDW is not the stable state of 7- 
Fe. 

The possibility of magnetic structures other than the 
1Q helical SDW has been suggested by Fujii et al. 28 on 
the basis of the LMTO calculations and the von Barth- 
Hedin potential. They compared the energies among 
the first-kind AF structure with Q — (1, 0, 0)27r/a, 
the commensurate 2Q structure with Q = (1,0, 0)27r/a 
and (0, 1, 0)27r/a, and the 3Q structure with Q = 
(l,0,0)27r/a, (0, 1, 0)2?r/a, and (0, 0, l)2?r/a. They con- 
cluded that the 3Q structure is the most stable at a = 
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6.8 a.u. (experimental lattice constant). Antropov and 
co- workers 29, 30 compared the energies of various mag- 
netic structures, including the structures obtained from 
the spin-dynamics calculations with 32 atoms in a unit 
cell. Using the GGA potential and the spdf basis, they 
claimed that the 3Q structure superposed with a heli- 
cal SDW with Q = (0,0,l/6)27r/o is the most stable 
for lattice constant a = 6.6 a.u., although it is nearly 
degenerate with a noncollinear eight-atom structure and 
helical structure. 

Kakehashi and coworkers 31 ' 32 have recently developed 
a molecular-dynamics (MD) approach which automati- 
cally determines the magnetic structure in a given unit 
cell at a finite temperature. Applying the approach to 
7-Fe with 500 atoms in a unit cell (5 x 5 x 5 fee lat- 
tice), 32 they found an incommensurate multiple spin den- 
sity wave (MSDW) state whose principal terms consist 
of 3Q MSDW with Q = (0.6, 0, 0)27r/a, (0, 0.6, 0)27r/a, 
and (0, 0, 0.6)27r/a. Subsequently, they performed the 
ground-state electronic-structure calculations 33 on the 
basis of the first-principles tight-binding LMTO method 
and the GGA potentials. Comparing the energies of var- 
ious magnetic structures, including the 1Q helical SDW 
and the MSDW found in the MD calculations, they 
concluded that the MSDW becomes the most stable 
state for lattice constants 6.8 < a < 7.0 a.u. More re- 
cently, Sjostedt and Nordstrom 27 implemented density- 
functional calculations based on the alternative lineariza- 
tion of the full-potential augmented-plane-wave method. 
They compared various competing collinear and non- 
collinear magnetic structures: ferromagnetic, commensu- 
rate 1Q, 2Q, and 3Q structures, and double-layered an- 
tiferromagnetic, as well as incommensurate helical struc- 
tures. For lattice constant a — 6.82 a.u., a collinear 
double-layered AF state was found to be the most stable 
structure. 

The results of the ground-state calculations described 
above indicate that we have not yet obtained a solid con- 
clusion as to the ground-state magnetic structure of cu- 
bic 7-Fe. In particular, no agreement between theory and 
experiment has yet been obtained. The discrepancy be- 
tween theory and experiment can be ascribed to both 
sides. On the experimental side, the difficulty originates 
in the fact that the bulk 7-Fe is stable only at high tem- 
peratures above the Curie temperature of a-Fe. At low 
temperatures, 7-Fe can be stabilized either in the form of 
small precipitates 5 in Cu or in thin Fe films 6, 15-19 grown 
on Cu. Thus, it is a subtle problem whether or not the 
SDW structure observed at low temperatures is identical 
with that should be realized in bulk cubic 7-Fe. Further- 
more, Tsunoda 5 emphasized in his detailed analyses that 
the 1Q helical SDW is one of possible structures with the 
same wave vector because of the high symmetry of the 
crystal structure of 7-Fe . 

On the theoretical side, the discrepancy should partly 
be ascribed to the various approximation schemes of the 
potential used in the density-functional calculations. As 
stated above, the equilibrium wave vectors for the IQ he- 
lical SDW predicted by the LDA calculations 20,21 are dif- 
ferent from those predicted by the GGA calculations. 22 
Similarly, the results of the atomic sphere approximation 



calculations 22 are not in agreement with those of full- 
potential calculations. 23-27 Through these calculations, 
it has become clear that there are various local minimum 
states in 7-Fe that are close in energy. 

In this situation, it is worthwhile to approach this 
problem from a phenomenological, but more general 
point of view so that one can gain useful information 
about possible magnetic structures of 7-Fe from only the 
symmetry of the system. The purpose of the present pa- 
per is to carry out such analysis by means of a Ginzburg- 
Landau type of theory and to clarify possible scenarios 
for the magnetic structure of 7-Fe from a phenomenologi- 
cal point of view. Since the observed magnetic moment 14 
of 7-Fe and calculated ones 28 near the equilibrium vol- 
ume are rather small (< 1/j.b), we will expand, in §2, a 
phenomenological free energy with respect to magnetic 
moments, and derive a general expression of free energy 
up to the fourth order. In §3 and 4, we discuss commensu- 
rate and incommensurate SDW structures, respectively, 
and obtain the magnetic phase diagrams in the space of 
expansion coefficients. The results of a preliminary anal- 
ysis for this part have been published. 34,35 In the last 
section, we discuss possible scenarios for the magnetic 
structure of 7-Fe that are consistent with the present 
theory. 

In contrast to the previous Landau type of phenomeno- 
logical theories applied to the incommensurate IQ SDWs 
and their harmonics in Cr 36 ' 37 and those applied to the 
commensurate MSDWs in 7-Mn alloys, 38 the present 
phenomenological free energy allows for incommensurate 
MSDWs with both linear and helical polarizations, in ad- 
dition to the commensurate MSDWs. In this respect, it 
should be emphasized that, so far, there has been no 
discussion about the possibility of helical MSDW states 
having the wave vectors found in the 7-FeCo precipitates 
in Cu, in spite of the fact that such MSDWs are also con- 
sistent with the experimental observation 5 of cubic 7-Fe. 

In the present analysis, we have shown that each 3Q 
state becomes the most stable state among the corre- 
sponding IQ, 2Q, and 3Q states for both commensurate 
and incommensurate wave vectors. This relation leads to 
magnetic phase diagrams that indicate the possibility of 
various 3Q MSDW states in fee transition metals, which 
are found to be consistent with the previous results: the 
commensurate 3Q state with Q = 2n/a for lattice con- 
stants 6.5 < a < 6.8 a.u. in the ground-state calcula- 
tions, 28,33 and the incommensurate linearly polarized 3Q 
state for a > 6.8 a.u. in the MD calculations. 33 The phase 
diagram also suggests a possibility of the incommensu- 
rate helically polarized 3Q state, which has not yet been 
investigated in the ground-state calculations. 

2. Phenomenological Free Energy 

We consider a free energy expansion with respect to 
the local magnetic moments on the fee lattice. Because 
the magnetic system in the absence of external magnetic 
fields has the time reversal symmetry, the free energy 
must include only even order terms with respect to mag- 
netic moments. The free energy per lattice site expanded 
up to the fourth order in magnetic moments can then be 
written as 
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1,1', l", V" a,0,f,S 

(2.1) 

Here, -/V is the number of lattice sites and to/ q denotes the 
a-component (a = x, y, z) of a magnetic moment on the 
Z-th site. a a p(l,l') and b a pjs(l,l' ,1" ,1"') are the expan- 
sion coefficients of the second- and fourth-order terms, 
respectively. Efc(En' I'M'") denotes the summations 
with respect to the site indices I and I' (I, I', I", and I'") 
over all integer values from 1 to N. Y^p" (EaS ^ 
denotes the summations with respect to the component 
indices a and /3 (a, f3, 7, and (5) over x, y, and z. 

Note that free energy (2.1) is written in the most gen- 
eral way in which each pair (quartet) of local magnetic 
moments at different sites is coupled via a coefficient 
independent of those for other pairs (quartets) of local 
magnetic moments. Therefore, free energy (2.1) can de- 
scribe the energy costs when the local magnetic moments 
change their magnitudes and directions, and hence can 
describe both the incommensurate and commensurate 
SDWs in itinerant magnets. The same type of free en- 
ergy expansion has been applied to the incommensurate 
SDWs of Cr. 36 < 37 

In the present paper, we apply free energy (2.1) to var- 
ious SDW states that can be realized on the fee transi- 
tion metals, specifically, 7-Fe. Then the free energy must 
be invariant with respect to the symmetry operations 
for the magnetic moments on the fee lattice: rotation 
TtiR, — ► 7£(m.R;), inversion m^ — ► m_H ; , and transla- 
tion rriRi — * mn l+r . Here, Ri is the position vector of 
the Z-th site; m^ = mi is the magnetic moment at the l- 
th site. 1Z denotes either the rotation C^lOO] or C3[lll], 
and r denotes an arbitrary lattice translation vector of 
the fee lattice. The requirement that the free energy be 
invariant under these operations yields the following free 
energy (see Appendix A): 



f = ^J2A(l,l')m l -m ll 



i,v 



+ E \B{l,l' ,1" ,l"'){mi ■ m v }{m vl ■ m vll } 

l,V,l",V" 

(y,z)(z,x)(x,y) 

+ C(l,l',l",l"') ]T (m, otYFll 1 aTTll" pTni" 't p 

+ mipm V f}rn v , a m l m a j\. (2.2) 

Here, A(l,l') = a xx (l,V), B(l,l',l",l"') = 
b xxxx (l, V, I", I'"), and C(l, V, I", V") = b yyzz (l, V , I", + 
b (1 1" V 1"') + b (I I" V" V\ - b (IV 1" 1"') 

u yzyz\ t 'j t ' 7 ^ ; ^ J ^ u yzzy\ l : 1 ; fc j 1 ) u xxxx\ t 'j'' 7 << 7 ^ )• 

The second-order terms with A(l,l') and the fourth- 
order terms with B(l,l',l" ,1"') in the free energy are 
isotropic since they are expressed in terms of the 



scalar products of magnetic moments. On the other 
hand, the fourth-order terms with the coefficients 
C(l,l',l",l"') are anisotropic. In the present paper, we 
restrict ourselves to the transition metals where the 
spin-orbit coupling effects are negligibly small, and thus 
neglect the anisotropic terms, i.e., we consider the case 
C(l, I', I", V") = in free energy (2.2). 

We now define the Fourier representation of the mag- 
netic moment at Ri by 



mi 



EBZ 

E 

<7 



m{q)e 



iq Hi 



(2.3) 



where X)q BZ denotes a summation with respect to q over 
the extended first Brillouin zone (EBZ) of the fee lattice, 
which is defined to include all the zone boundary points. 
This form, with the use of the EBZ, has the merit that 
one can argue the structures in both the commensurate 
and incommensurate cases on the same footing. 

The Fourier representation of the isotropic free energy 
is then given by 



EBZ 

f = J2 A(Q)\m(q)\ 2 + ]T E B (l> q", O 

q K q,q',q",q"' 

x {m(q) ■ m(q')}{m(q") ■ m(q"')}, (2.4) 

where the coefficients A(q) and B(q, q' , q" ,q"') are de- 
fined by 



(2.5) 



B(q,q',q",q"') = ^s E B(n,n',n") 



xe iq'-R ne i q "-R n , e iq'"-R n „ S 



q+q'+q"+q'",K- 
(2.6) 

In defining the coefficients A(q) and B(q, q' , q" , q'") in 
cqs. (2.5) and (2.6), we have introduced the relative 
coordinates R n = Rv — Ri, R n ' = Ri» — Ri, and 
R n " = Ri 1 " — Ri, and have used the notations 

A(n) = A(l, I') = A{Ri,Rp) = A(0, R n ), (2.7) 

B(n, n', n") = B(l, I', I", I'") = B(R U R V , R v , , R v „ ) 

= B(0, R n , R n i,R n „). 

(2.8) 

Here, the last equalities in eqs. (2.7) and (2.8) result 
from the translational symmetry of the fee lattice. Note 
that the inclusion of the fourth-order terms in free energy 
(2.1), and hence in eq. (2.4), is essential for the present 
analysis since the second-order terms alone do not de- 
scribe the MSDW states. 

As is always the case with all Landau- type theories, 
one should keep in mind the applicability of free energy 
expansion (2.4), which is valid for the system with small 
local magnetic moments. In the present paper, our main 
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concern is the complex magnetic structures of 7-Fe which 
appear between the nonmagnetic state and the strong 
ferromagnetic state with increasing volume. In such a 
region, the magnitudes of the magnetic moments are rel- 
atively small. Even if this were not the case, we can apply 
the theory at any volume near the transition tempera- 
ture where the local magnetic moments become small. 
Because of these reasons, we apply free energy expan- 
sion (2.4) for the analysis of various SDW states in 7-Fe, 
i.e., the commensurate 1Q, 2Q, and 3Q SDWs, and the 
incommensurate linearly and helically polarized SDWs 
with 1Q, 2Q, and 3Q wave vectors. On the basis of the 
analysis, one can draw a general and exact conclusion 
independent of the specific model or the potential in the 
first-principles calculations. 

3. Commensurate SDW Structures 

We investigate here the commensurate SDW struc- 
tures whose magnetic moments are given by 

m t = J2[m(Q n )e^- R ' + m(Q n )e-^- R <], (3.9) 

n=l 

with the set of equivalent wave vectors Q l = 
(l,0,0)(27r/a), Q 2 = (0, 1, 0)(27r/o), and Q 3 = 
(0, 0, l)(2n/a). Here, m(Q 1 ), m(Q 2 ), and m(Q 3 ) are 
real and assumed to be orthogonal to each other: m(Q 2 ) • 
m(Q 3 ) = m(Q 3 ) ■ m(Qi) = m{Q x ) • m(Q 2 ) = 0. The 
commensurate MSDW with the form of eq. (3.9) has 
been discussed in the previous ground-state electronic- 
structure calculations. 27-29,33 The free energy is given 

by 



fco = ]T[i Q |m(Q 4 )| 2 + (B 1Q + B 2 Q)\m(QA\ 4 } 

i=i 

(2,3)(3,1)(1,2) 

+ £ Bxoolm^lV^OI 2 - (3.10) 

The coefficients Aq, B\q, B 2 q, and B\qq are expressed 
in terms of linear combinations of the coefficients {^4(q)} 
and {B(q, q', q", q'")} in cqs. (2.5) and (2.6), with q, q' , 
q", and q'" being chosen from ±Q 1 , ±Q 2 , and ±Q3- 
The full expressions of these coefficients are given in Ap- 
pendix B. We see that free energy (3.10) depends only 
on the absolute squares of magnetic moments |m(Qj)| 
(i = 1,2,3). Thus the SDW structures described by 
this free energy are degenerate with respect to the di- 
rections of polarization. This degeneracy is partially re- 
moved when the anisotropic terms are included in the 
free energy. In the following, we investigate three com- 
mensurate magnetic structures, the first-kind AF struc- 
ture, the 2Q structure, and the 3Q structure, on the 
basis of free energy (3.10). 

3.1 First-kind antijerromagnetic structure 

The equilibrium magnetic moment for the first-kind 
AF structure is obtained by minimizing free energy 
(3.10) with respect to \m{Q 1 )\ 2 and setting m(Q 2 ) = 



m(Q 3 ) = 0. We have 
under the condition 



Ac 



1 1/2 



2{B 1Q + B 2Q ) 



A 



Q 



> 0. 



(3.11) 



(3.12) 



2{B lQ + B 2 q) 

In order for the solution (3.11) to be thermodynamically 
stable, it is necessary that 



d 2 f 



> 0. 



(3.13) 



(3.11) 



Inequalities (3.12) and (3.13) reduce to 

A Q < 0, (3.14) 

B 1Q +B 2Q >{). (3.15) 

Inequalities (3.14) and (3.15) yield the stability condi- 
tions for the first-kind AF structure. 

The equilibrium free energy is obtained by substituting 
eq. (3.11) and m(Q 2 ) = m(Q 3 ) = into eq. (3.10): 



/af 



4 



(3.16) 



4(5iq + B 2Q ) 

The amplitude M of the magnetic moment per site is 
given by 



N 

EBZ 



mi ■ mi 



(3.17) 



q,q> K 

In the commensurate case, it becomes 



M 



4(|m(Q 1 )| 2 + |m(Q 2 



|m(Q 3 )| 2 ). (3.18) 



Substituting eq. (3.11) and m(Q 2 ) = m(Q 3 ) = into 
eq. (3.18), we have the amplitude Maf of the magnetic 
moment for the first-kind AF structure: 



M 



2A, 



AF 



Q 



Biq + B 2Q 



(3.19) 



3.2 2Q structure 

The equilibrium magnetic moments for the 2Q struc- 
ture are obtained by minimizing free energy (3.10) with 
respect to |m.(<5i) | 2 an d |m(Q 2 )| 2 and setting m(Q 3 ) = 
0, we have 

A Q + 2(B 1Q + B 2 g)|m(0i)| 2 + B 1QQ \m(Q 2 )\ 2 = 0, 



(3.20) 

A Q + 2(B 1Q + B 2Q )\m(Q 2 )\ 2 + B^miQ^ 2 = 0. 



(3.21) 
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When 

D 2Q = 4 (^iQ + B2Q? - B\ QQ + 0, 
eqs. (3.20) and (3.21) are solved as 

= |m(Q 2 )| = 
under the condition 



(3.22) 



1/2 



2 (- B iQ + B 2 q) + Biqq 



(3.23) 



A 



Q 



> 0. 



(3.24) 



2(-Biq + B 2Q ) + Biqq 

In order for the solution (3.23) to be thermodynamically 
stable, it is necessary that 



2 2 



* 2 / = EE 



-ij-idUmiQ^jdUmiQ^} 



(3.23) 



x ( 5|m(g j )| 2 5|m(Q j )| 2 >0. (3.25) 
This condition is equivalent to 



fn > 0, 
where fa is defined by 



fn fi2 
/21 /22 



>0, 



fij 



d 2 f 



dUmiQ^jdUmiQ^} 



(3.23) 



Using eq. (3.27), condition (3.26) becomes 
2(B 1Q + B 2Q ) >0, 

\2 52 



(3.26) 

= 1,2). 
(3.27) 



(3.28) 

4(Si + B 2Q f - B{qq > 0. (3.29) 
Conditions (3.22), (3.24), (3.28), and (3.29) reduce to 

A Q < 0, (3.30) 
\Biqq\ 



B\Q + B 2 Q > 



(3.31) 



Inequalities (3.30) and (3.31) yield the stability condition 
for the 2Q structure. 

The equilibrium free energy is obtained by substituting 
cq. (3.23) and m(Q 3 ) = into eq. (3.10): 



A 2 



2{B 



IQ 



B 



2QJ 



Biqq 



(3.32) 



The amplitude M 2 q of the magnetic moment per site is 
obtained by substituting eq. (3.23) into eq. (3.18): 



M 2 - 

2Q 



$Aq 



2(B 1Q + B 2Q ) + B 1QQ 



(3.33) 



3.3 3Q structure 

The equilibrium magnetic moments for the 3Q struc- 
ture are obtained by minimizing free energy (3.10) with 
respect to ^(Q^] 2 , |m(Q 2 )| 2 , and |m(Q 3 )| 2 . We ob- 
tain 



MQJI = \m(Q 2 )\ = \m(Q 3 )\ 

Aq 



2(B 1Q + B 2Q + Biqq) 



1/2 



under the condition 



Ac 



> 0. 



(3.34) 



(3.35) 



2(B 1Q + B 2Q + B 1QQ ) 

The thermodynamical stability analysis and inequality 
(3.35) lead to the stability condition for the 3Q structure: 



An <0, 



B\q + B 2 q > 



B 



1QQ 



(3.36) 

for Biqq > 0, (3.37) 



Biq + B 2 q > —Biqq for B 1QQ < 0. (3.38) 

The equilibrium free energy is obtained by substituting 
cq. (3.34) into eq. (3.10): 

Uq 1 



3^ 



• (3.39) 

4(Siq + B 2Q + Biqq) 

The amplitude M 3 q of the magnetic moment for the 

3Q structure is obtained by substituting eq. (3.34) into 
cq. (3.18): 

GAq 



M A ~ = -- 

3Q 



B\Q + B 2 q + Biqq 



(3.40) 



3.4 Relative stability among commensurate structures 

The relative stability among the first-kind AF, and the 
2Q and 3Q structures has been determined by compar- 
ing stability conditions (3.14), (3.15), (3.30), (3.31), and 
(3.36)-(3.38), and equilibrium free energies (3.16), (3.32), 
and (3.39). The obtained magnetic phase diagram 39 for 
Aq < is shown in Fig. 1 in the space of expansion co- 
efficients Bxqq/Biq and B 2Q /Bi Ql where B Y q > for 
B 2Q /B 1Q > -1 and B 1Q < for B 2Q /B 1Q < -1. 

All three commensurate structures appear in the mag- 
netic phase diagram. In the AF phase (0 < Biq + B 2 q < 
I Biqq |/2), the AF state is the only stable structure. In 
the 2Q phase (0 < -B 1QQ /2 < B 1Q + B 2Q < -B 1QQ ), 
the AF and 2Q structures are stable. Comparison of 
their free energies shows that the 2Q structure is the 
most stable state in this region. Note that the ampli- 
tude of the magnetic moment for the 2Q structure is 
larger than that for the AF structure in this region. In 
the 3Q phase (0 < B 1QQ /2 < B 1Q + B 2Q , < -B 1QQ < 
B\q + B 2 q), all three commensurate structures are sta- 
ble. Since the equilibrium free energies satisfy the in- 
equality /af > f 2 Q > / 3 q, the 3Q structure is the most 
stable state in this region. The relation among the ampli- 
tudes of the magnetic moments for the three structures 
Maf < M 2 n < M„q shows that the 3Q structure is the 



state with the largest magnetic moment in this region. 

It is of interest here to compare the present magnetic 
phase diagram for the commensurate structures with 
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Fig. 1. Magnetic phase diagram for the commensurate structures 
with Q = 27r/a for Aq < 0. The first-kind antiferromagnetic 
(AF), and the 2Q, and 3Q phases are shown in the space of 
expansion coefficients Biqq / 'B\q and B 2 q / B\q, where B\q > 
for B 2Q /B 1Q > -1 and S 1Q < for B 2Q /B 1Q < -1. 

the past results of the ground-state electronic-structure 
calculations for cubic 7-Fe. Those calculations 20-27 ' 32,33 
showed that the magnetism of 7-Fe depends sensitively 
on the volume; the first-kind AF state appears for lat- 
tice constants a < 6.5 a.u., complex magnetic structures 
for lattice constants 6.5 < a < 7.0 a.u., and the fer- 
romagnetic state for 7.0 a.u. < a. Although there is a 
wide diversity in the results of the predicted magnetic 
structures for intermediate values of the lattice constant, 
6-5 < a < 7.0 a.u., it is worth noting that the ground- 
state electronic-structure calculations by Kakehashi et 
al. 33 and those by Fujii et al. 28 predicted the commensu- 
rate 3Q state to appear for lattice constants a < 6.8 a.u. 
Kakehashi et al. 33 performed the first-principles tight- 
binding LMTO calculations for 7-Fe using the GGA po- 
tential and compared the ground-state energies of various 
MSDW states. It was predicted that with increasing vol- 
ume, 7-Fe undergoes a transition from the first-kind AF 
state to the commensurate 3Q state at a = 6.5 a.u. with 
the 3Q state remaining stable until a — 6.8 a.u. This 
result is consistent with the present magnetic phase dia- 
gram of Fig. 1 which shows the possibility of the transi- 
tion from the AF to the 3Q phase crossing the boundary. 

Fujii and co-workers 28 found three possible ground 
states for 7-Fe: the first-kind AF, and the 2Q and 3Q 
structures, using the LMTO method and the von Barth- 
Hedin potential, and found numerically that the 3Q 
structure is the most stable among the three at a = 6.8 
a.u. As we mentioned above, the present theory yields 
the same relative stability, /af > J 2 q > /3Q1 among 
the AF, 2Q, and 3Q structures; their results are veri- 
fied by the present theory. It is of interest to note that 
the amplitudes of their magnetic moments for the 3Q 



and 2Q structures were found to be the same and larger 
than that for the AF structure: M 3 q = M 2 q > Maf- 
According to cqs. (3.33) and (3.40), this implies that 
B\q + B 2 q = Biqq/2; the ground state of 7-Fe cal- 
culated by Fujii et al. is located in the vicinity of the 
AF-3Q boundary in the 3Q phase in Fig. 1. 

4. Incommensurate SDW Structures 

We consider SDW structures described by three in- 
commensurate wave vectors Q lt Q 2 , and Q 3 . These wave 
vectors are assumed to be equivalent in space with each 
other and to satisfy the following incommensurate con- 
ditions: 

±4Q^K, ±2Qi?K (2 = 1,2,3), (4.41) 

±2(Q i ±Q j )^K, ±(3Q i ±Q j )^K, 

±(Q l ±Q J )^K ((i,j) = (2,3)(3,l)(l,2)), (4.42) 

± (2Q, ± Qj ±Q k )^K 

((i,j,fc) = (l,2,3)(2,3,l)(3,l,2)). (4.43) 

Here, K is a reciprocal lattice vector of the fee lattice. 
Incommensurate conditions (4.41)-(4.43) are satisfied by 
the wave vectors predicted for the bulk cubic 7-Fe in 
the electronic band-structure calculations 20 ' 21 and in the 
molecular-dynamics calculations, 32 ' 33 and by that found 
in the neutron diffraction measurements 5 of the cubic 7- 
Feioo-ajCoa; (x < 4) precipitates in Cu. The free energy 
describing the incommensurate SDWs can be written as 

3 

/ ic = ^[yl Q |m(QJ| 2 + J B 1Q |m(Q i )| 4 

i=i 

+ B 2Q m 2 (Q i )m* 2 (Q i )] 

(2,3)(3,1)(1,2) 

+ Yl [BiQQ\m(Q i )\ 2 \m(Q j )\ 2 

+ B 2QQ \m(Q l )-m(Q J )\ 2 + B 3QQ \m(Q l )-m*(Q 1 )\ 2 ]. 

(4.44) 

The coefficients Aq, B 1Q , B 2Q , B 1QQ , B 2QQ , and B 3QQ 
are expressed in terms of linear combinations of coeffi- 
cients A(q) and B(q,q',q",q"') in cqs. (2.5) and (2.6), 
with q, q', q" , and q'" being chosen from Q 1 , Q 2 , and 
Q 3 satisfying conditions (4.41)-(4.43). The full expres- 
sions of these coefficients are given in Appendix B. On 
the basis of free energy (4.44), we investigate the 1Q, 
2Q, and 3Q linearly polarized SDWs, and the 1Q, 2Q, 
and 3Q helically polarized SDWs. 

4-1 Linearly polarized SDWs 

We consider first the linearly polarized SDWs whose 
magnetic moments are described by 

3 

mi = ^[m(QJe i£ ^ + m*{Q n )e-^ R % (4.45) 

n=l 
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with 

m(QJ = (m x (Q n ),m y (Q n ),m z (Q n ))e ia " (n = 1,2,3). 

(4.46) 

Here, m x (Q n ), m y (Q n ), and m z {Q n ) (n = 1,2,3) are 
assumed to be real. ati, a 2 , and «3 are phase factors. We 
consider the case in which m(Q 1 ), m(Q 2 ), and m(Q 3 ) 
are orthogonal to each other: 

m(Q 2 )-m(Q 3 ) = m(Q 3 )-m(Q 1 ) = m(Q 1 )-m(Q 2 ) = 0. 

(4.47) 

The free energy for the linear SDWs is obtained by 
substituting eqs. (4.46) and (4.47) into eq. (4.44): 



3 

E 



[A Q \m(Q t )\ 2 + (B 1Q + B 2Q )\m(Q t )\ 4 



(2,3)(3,1)(1,2) 

+ J2 S 1QQ |m(Q i 



lV(Q,)l 2 - 



(4.48) 



Note that free energy (4.48) depends only on the abso- 
lute squares of magnetic moments, \m(Q 1 )\ 2 , |m(Q 2 )| 2 , 
and |m(Q 3 )| 2 . This again implies that the IQ, 2Q, and 
3Q states are degenerate with respect to the directions of 
polarization. This degeneracy is partially removed when 
the anisotropic terms are included in the free energy. We 
also note that free energy (4.48) has the same form as 
eq. (3.10) in which Aq, B 2 q, and B\qq have been re- 
placed by Aq, B 2 q, and Biqq, respectively. Therefore, 
following the same steps as in §3, we obtain the equi- 
librium states of the IQ, 2Q, and 3Q linear SDWs as 
follows. 

4-1.1 1Q linearly polarized SDW 
Magnetic moment 



MQi)l = 

Stability condition 



Ac 



^{BiQ + B 2 q) 



1/2 



Aq < 0, 
B\q + B 2 q > 0. 
Equilibrium free energy 



fiQ = ~- 



A 2 



4(Si Q + B 2Q ) ' 
Amplitude of the magnetic moment 



M\ Q = 2\m(Q 1 ) 



A ( 



BlQ + B 2 Q 



4-1-2 2Q linearly polarized SDW 
Magnetic moment 



\m(Q,)\ = \m(Q 2 )\ = 
Stability condition 



Ac 



2(B 1Q + B 2Q )+B 



iQQ 



Aq < 0, 



(4.49) 

(4.50) 
(4.51) 

(4.52) 
(4.53) 

- 1/2 

(4.54) 
(4.55) 



B\q + B 2 q > 



\Biqq\ 



Equilibrium free energy 



2 (-BiQ + B 2 q) + Biqq 
Amplitude of the magnetic moment 

M\ 4Aq 

2Q 2 (-BlQ + B 2Q ) + Biqq 

4-1-3 3Q linearly polarized SDW 
Magnetic moment 



(4.56) 



(4.57) 



(4.58) 



\m{Q x )\ = \m(Q 2 ) 



\m(Q 3 )\ 
Aq 



2{Biq + B 2 q + Biqq) 



1/2 



Stability condition 
Aq < 0, 

d , r ^ Biqq 
o\q + o 2 r> > — - — 



B\Q + B 2 q > -Biqq 

Equilibrium free energy 
/3Q 



for B 1QQ > 0, 
for B 1QQ < 0. 



3^> 



4(Biq + B 2 q + B 1QQ ) 
Amplitude of the magnetic moment 



3A q 



BlQ + B 2Q + Biqq 



(4.59) 

(4.60) 
(4.61) 
(4.62) 

(4.63) 
(4.64) 



4-2 Relative stability among linear SDWs 

The relative stability among the incommensurate 1Q, 
2Q, and 3Q linear SDWs has been determined by com- 
paring stability conditions (4.50)-(4.51), (4.55)-(4.56) 
and (4.60)-(4.62), and equilibrium free energies (4.52), 
(4.57), and (4.63). The obtained magnetic phase diagram 
for Aq < is shown in Fig. 2 in the space of expansion 
coefficients Biqq/Biq and B 2Q /Bi Q , where B 1Q > 
for B 2Q /B 1Q > -1 and B 1Q < for B 2Q /B 1Q < -1. 

The relative stability among the linear SDWs has been 
found to have the same feature as that of the rela- 
tive stability among the commensurate SDWs. In the 
IQ phase (0 < B 1Q + B 2Q < \B 1QQ \/2), the 1Q lin- 
ear SDW is the only stable structure. In the 2Q phase 
(0 < -B 1QQ /2 < Biq + B 2Q < —Biqq), both the IQ 
and 2Q linear SDWs are stable, but the latter has a lower 
free energy and larger amplitude of the magnetic mo- 
ment. In the 3Q phase (0 < B 1QQ /2 < B 1Q + B 2Q , < 
~Biqq < Biq + B 2 q), all three linear SDWs are sta- 
ble, and the 3Q state has the lowest free energy and the 
largest amplitude of the magnetic moment. 

Concerning the magnetism of 7-Fe, the present mag- 
netic phase diagram for the linear SDWs indicates the 
possibility of the 3Q and 2Q MSDWs as well as the 
IQ SDW. In the ground-state calculations of bulk cu- 
bic 7-Fe, although the magnetic structure for lattice 
constants 6.5 < a < 7.0 a.u. is under debate, the 
possibility of incommensurate linear MSDW was sug- 
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Fig. 2. Magnetic phase diagram for the incommensurate linear 
SDWs for A Q < 0. The 1Q, the 2Q, and the 3Q phase are 
shown in the space of expansion coefficients Biqq/Biq and 
B 2 q/B 1q , where B 1Q > for B 2Q /B 1Q > -1 and B 1Q < 
for B 2Q /B 1Q < -1 . 



Here, (j,k,m) is (1,2,3), (2,3,1), and (3,1,2) when j = 
1, 2, and 3, respectively, ei, e 2 , and e 3 form an ortho-nor- 
mal basis set. a\, a 2 , and 03 are phase factors. Note 
that the possibility of the helical 3Q MSDW given by 
cq. (4.65) has not been examined in either the experi- 
mental analyses or the electronic-structure calculations. 

Introducing the basis set (ejk, ej m ) obtained by a rota- 
tionof(e fc ,e m )bya, {{j,k,m) = (1, 2, 3)(2, 3, 1)(3, 1, 2)) 
for each helical component j, 

&jk = efc cos aj + e m sin aj , (4.66) 

6j m — sin aj + e m cos aj , (4.67) 

we have the expression 

3 

mi = Y^[m(Qj)e iQ r R ' + m*(Qj)e- iQ r R '}. (4.68) 



Here, 



m(Q,) = 



= (1,2,3)(2,3,1)(3,1,2)). (4.69) 



The free energy for the helical SDWs is obtained by 
substituting eq. (4.69) into eq. (4.44): 



gested by Kakehashi and coworkers 32,33 On the basis 
of the molecular-dynamics (MD) method, 32 they pre- 
dicted a new MSDW whose principal terms consist of 
3Q waves with Q = (0.6, 0, 0)(27r/a), (0, 0.6, 0)(27r/a), 
and (0, 0, 0.6)(27r/a). Subsequently, they performed the 
ground-state electronic-structure calculations 33 using the 
first-principles tight-binding LMTO method and the 
GGA potentials to compare the ground-state energies 
of various magnetic structures: the first-kind AF state, 
the commensurate 3Q structure, the incommensurate 1Q 
helical SDW, the incommensurate MSDW found in the 
MD calculations, and the ferromagnetic state. It was con- 
cluded that the MSDW becomes the most stable state 
for lattice constants 6.8 < a < 7.0 a.u. In particular, 
they find that the incommensurate 3Q MSDW is sta- 
ble as compared with the 1Q SDW irrespective of the 
lattice constant and that the amplitude of the magnetic 
moment for the 3Q state is larger than that for the IQ 
state. These results are consistent with the present re- 
sult that the 3Q MSDW is always stabilized and has a 
larger amplitude of the magnetic moment as compared 
with the IQ SDW when the 3Q solution exists, resulting 
in a wide range of the 3Q phase in the magnetic phase 
diagram in Fig. 2. 

4-3 Helically polarized SDWs 

Next, we consider the helically polarized SDWs whose 
magnetic moments are described by 



mi = ^2 V2\m(Qj)\[e k cos(Qj ■ Ri + aj) 

3=1 



+ e m sm(Q r Ri+a 3 )}. (4.65) 



f K = J2[A Q \m(Q l )\ 2 + B 1Q \m(Q l )\ 4 



i=l 



(2,3)(3,1)(1,2) 

+ E (Bi QQ +B 2QQH )\m(Q l )\ 2 \m(Q J 

(i,3) 



with 



B 



2QQH 



_ B 2 QQ + B 3 QQ 



(4.70) 



(4.71) 



Note that free energy (4.70) depends only on the absolute 
squares of magnetic moments |m(Q 1 )| 2 , |m(Q 2 )| 2 , and 
\m(Q 3 )\ 2 ; therefore, the IQ, 2Q, and 3Q helical states 
are degenerate with respect to the directions of polariza- 
tion. We also note that free energy (4.70) is identical to 
cq. (3.10) in which Aq, B\q + B 2 q, and Biqq have been 
replaced by A Q , B X q, and Biqq + B 2 qqk- Thus, follow- 
ing the same steps as in §3, we obtain the equilibrium 
states of IQ, 2Q, and 3Q helical SDWs as follows. 

4-3.1 IQ helically polarized SDW 
Magnetic moment 



MQi)l - 

Stability condition 



A 



Q 



2B 



iQ 



1/2 



Aq <0, 
Biq > 0. 



(4.72) 



(4.73) 
(4.74) 



Equilibrium free energy 
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/lQH 



Amplitude of the magnetic moment 



1Q 

4.3.2 2Q helically polarized SDW 
Magnetic moment 



MQOI - \m(Q 2 )\ = 
Stability condition 



(4.75) 



(4.76) 



Aq 



2B 1Q + B 1QQ + B 2 qqh 



Aq<0, 

Biq > -^\B\qq + B 2 qqh\. 



Equilibrium free energy 
/2QH = - 



A% 



2Bi Q + Biqq + B 2 qqu 
Amplitude of the magnetic moment 



4A 



Q 



2B 1Q + Biqq + B 2 qqh 



1/2 

(4.77) 

(4.78) 
(4.79) 

(4.80) 

(4.81) 



4.3.3 3Q helically polarized SDW 
Magnetic moment 



MQi)! - \m(Q 2 )\ = |ro(Q 3 )| 



2 (Biq + Biqq + B 2 qq h ) 

Stability condition 
Aq <0, 

Biqq + B 2Q qh 



1/2 



(4.82) 



(4.83) 



Biq > 



for Biqq + B 2 q QH > 0, 
(4.84) 



Biq > -{B 1QQ + B 2QQ k) for B 1QQ + B 2 qqu <{). 

(4.85) 

Equilibrium free energy 
/3QH 



*A% 



4(-Biq + Biqq + B 2 qqh) 
Amplitude of the magnetic moment 



M| QH = J Aq - 

Biq + Biqq + B 2QQ n 



(4.86) 



(4.87) 




Fig. 3. Magnetic phase diagram for the helical SDWs for Aq < 
and Biq > 0. The 1Q helical (1QH), 2Q helical (2QH), and 
3Q helical (3QH) phases are shown in the space of expansion 
coefficients B 1Q q/B 1q and B 2 qqk/B 1 q. 



4-4 Relative stability among helical SDWs 

The relative stability among the incommensurate 1Q, 
2Q, and 3Q helical SDWs has been determined by com- 
paring stability conditions (4.73)-(4.74), (4.78)-(4.79) 
and (4.83)-(4.85), and the equilibrium free energies 
(4.75), (4.80), and (4.86). The obtained magnetic phase 
diagram is shown in Fig. 3 in the space of expansion co- 
efficients Biqq/Biq and B 2 qq^/B 1q for A Q < and 
Biq > for which the solutions of helical SDWs exist. 

The relative stability among the helical SDWs has 
been found to have the same feature as that of the 
relative stability among the commensurate SDWs and 
among the linear SDWs. In the 1Q phase (0 < Biq < 
\B\qq + B 2QQ h|/2), the 1Q helical SDW is the only 
stable structure. In the 2Q phase (0 < ~{B\qq + 
B 2 qqh)/2 < Biq < -(Biqq + B 2QQH )), both the 1Q 
and 2Q helical SDWs are stable, but the latter has a 
lower free energy and larger amplitude of the magnetic 
moment. In the 3Q phase (0 < (Biqq + B 2 qqh)/2 < 
Biq, < -(Biqq + B 2 qqh) < Biq), all the three heli- 
cal SDWs arc stable, but the 3Q state yields the lowest 
free energy and the largest amplitude of the magnetic 
moment. 

Neutron diffraction experiments 5 on cubic 7- 
Feioo-ajCOa; (a; < 4) alloy precipitates in Cu 
showed a magnetic satellite peak for wave vector 
Q = (0.1, 0, l)2n/a. The magnetic structure was sug- 
gested to be a helical SDW but has not been determined 
precisely. This is because the neutron diffraction analysis 
cannot distinguish between the 1Q and 3Q states 40 
when the crystal structure of the 7-Fe precipitates 
is properly cubic and the distribution of domains is 
isotropic. The present finding that the 3Q helical 
MSDW is always stable as compared with the 1Q and 
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2Q SDWs when the 3Q solution exists suggests that the 
3Q helical MSDW should be taken into consideration in 
addition to the 1Q helical SDW in the analysis of the 
magnetic structure of cubic 7-Fe. 

One might have a question as to why the 3Q MSDW 
is stable in a wide range of the parameter region (4.84) 
or (4.85), while the previous phenomenological theories 
concerning the Heisenberg model 41-43 predict the 1Q he- 
lical SDW ground state. The physical reason for this is as 
follows. For simplicity, we consider first the free energy 
for the helical SDW states eq. (4.70) without the mode- 
mode coupling term (the term with {B\qq + B^qqh))- 
The free energy for the 3Q state is three times smaller 
than that for the 1Q state, as is seen from eqs. (4.75) and 
(4.86). This free energy gain is caused by the increase in 
amplitudes of local magnetic moments as is seen from 
eqs. (4.76) and (4.87). This is characteristic of the itin- 
erant electron system. In the localized model reported 
by Yoshimori 41 and Kaplan, 42,43 this mechanism of en- 
ergy gain is forbidden because of the constraint of the 
constant amplitudes of local magnetic moments, so that 
the 1Q state is realized. Under the constraint of a con- 
stant amplitude of local magnetic moments, the present 
theory also predicts the 1Q helical SDW as the stable 
structure, which is consistent with the theory presented 
by Yoshimori and Kaplan. 

When the mode-mode coupling term is positive, it sup- 
presses the increase in the amplitudes of local moments 
of the 3Q state (see eqs. (4.76) and (4.87)). As a result, 
the 3Q MSDW is stable only when the coefficient of the 
mode-mode coupling term is smaller than a critical value; 
otherwise, the 1Q helical SDW is stable. This condition 
is given by inequality (4.84). 

4-5 Relative stability among linear and helical SDWs 

In the previous two subsections, we examined the rel- 
ative stability among the linear SDWs and that among 
the helical SDWs, separately, assuming incommensurate 
conditions (4.41)-(4.43) for the wave vectors. Free en- 
ergy (4.44) having such incommensurate wave vectors, 
however, allows for both linear and helical SDWs in the 
common space of the expansion coefficients. In order to 
discuss their relative stability, we present, in this subsec- 
tion, a magnetic phase diagram allowing for both linear 
and helical SDWs. 

Comparing the equilibrium free energies for the linear 
and helical SDWs, we have obtained magnetic phase dia- 
grams 44 for Aq < and Biq > for which both the lin- 
ear and helical SDWs are stable. Figure 4 shows an exam- 
ple of the magnetic phase diagram for B^qq^I B\q = 1 
in the space of expansion coefficients Biqq/Biq and 
B2Q/B1Q, where Aq < and Biq > 0. We see that the 
3Q linear (3Q) and 3Q helical (3QH) MSDWs occupy 
most of the region —2 < Biqq/Biq < 1. This arises from 
the fact that the 3Q state is stable when the mode-mode 
coupling term Biqq or Biqq + B 2 qqh is relatively small, 
as discussed in §4.4. Although we presented an example 
of the magnetic phase diagram for B 2 qqy{I B\q = 1 in 
Fig. 4, changing the value of B 2 qqh/Biq does not alter 
the global feature of the magnetic phase diagram if it 
does not become exceedingly large. Thus we discuss the 
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Fig. 4. Magnetic phase diagram for the incommensurate SDWs 
for Aq < 0, Biq > and Bzqqyi/ Biq = 1. The phases of the 
1Q linear SDW (1Q), the 2Q linear MSDW (2Q), the 3Q linear 
MSDW (3Q), the 1Q helical SDW (1QH), the 2Q helical MSDW 
(2QH), and the 3Q helical MSDW (3QH) arc shown in the space 
of expansion coefficients Biqq/Biq and B2Q/B1Q. Coexistence 
lines between the linear and helical SDWs are indicated by solid 
lines. 

possible magnetic structures of 7-Fe on the basis of the 
magnetic phase diagram in Fig. 4. 

According to the ground-state electronic-structure cal- 
culations by Kakehashi et al. 33 the 3Q linear MSDW is 
stabilized for lattice constants 6.8 < a < 7.0 a.u. This 
result can be explained by the existence of the 3Q lin- 
ear MSDW phase in Fig. 4. Note that the 3Q linear 
MSDW solution is extended to the region of the 3Q he- 
lical MSDW in Fig. 4. Thus there is another possibil- 
ity that the latter MSDW is stable when the 3Q helical 
MSDW is taken into account in the electronic-structure 
calculations. It is highly desirable to investigate the rel- 
ative stability between 3Q linear and helical MSDWs in 
the first-principles ground-state calculations of 7-Fe. 

5. Summary and Discussion 

We have investigated the relative stability among var- 
ious SDW structures in fee transition metals on the basis 
of a Ginzburg-Landau type of free energy with the terms 
up to the fourth order in magnetic moments. We have 
obtained magnetic phase diagrams in the space of ex- 
pansion coefficients for both commensurate and incom- 
mensurate wave vectors, and discussed their implications 
on the magnetism of cubic 7-Fe. 

In both the commensurate and incommensurate cases, 
we proved that the 3Q state is always stable compared 
with the corresponding 2Q and 1Q states when there is a 
solution of the 3Q state. The energy gain of the 3Q state 
is caused by a change in the amplitude of the magnetic 
moment when an additional mode Q is introduced. Ac- 
cordingly, we have the relation M 3 q > M 2 q > M\q. This 
is characteristic of itinerant magnets and has a profound 
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effect on the magnetic phase diagram. In the localized 
systems, only the 1Q helical SDW is possible because of 
the fixed amplitudes of local magnetic moments. 

On the basis of the magnetic phase diagrams for the 
commensurate case (Fig. 1) and for the incommensu- 
rate case (Figs. 2-4), we have discussed the possible 
magnetic structures of cubic 7-Fe. According to the 
ground-state electronic-structure calculations, 20-27 ' 32 ' 33 
the magnetism of 7-Fe depends sensitively on the vol- 
ume; the first-kind AF state appears for lattice constants 
a < 6.5 a.u., SDW structures for lattice constants 6.5 
< a < 7.0 a.u., and the ferromagnetic state for a > 7.0 
a.u. The magnetic structures for 6.5 < a < 7.0 a.u. are 
under debate and there is a wide diversity in the pre- 
dicted results. 

The ground-state electronic-structure calculations by 
Kakehashi et al. 33 and those by Fujii et al. 28 predicted 
the commensurate 3Q state to appear for lattice con- 
stants a < 6.8 a.u. This result can be explained by the 
existence of the commensurate 3Q phase, as shown in 
Fig. 1, specifically in the region of the phase diagram 
with B 1QQ w 2(B 1Q + B 2Q ) and B 1QQ > 0. 

For larger lattice constants, 6.8 < a < 7.0, the 
MD calculations 33 for 7-Fe predicted the incommensu- 
rate 3Q linear MSDW state with Q = (0.6, 0, 0)27r/a, 
(0, 0.6, 0)27r/a, and (0, 0, 0.6)27r/a. This can be explained 
by the existence of the 3Q linear MSDW phase in Fig. 4. 
It should be noted, however, that there is another possi- 
bility of the 3Q helical phase since the 3Q linear MSDW 
solution is extended to the region of the 3Q helical phase. 
It is desirable to examine the relative stability between 
the 3Q linear and helical states for the above wave vector 
by means of the ground-state electronic-structure calcu- 
lations. 

Experimentally, the SDW of cubic 7-Fe was found 
for wave vector Q = (0.1, 0, l)27r/a. 5 It was suggested 
that a helical spin configuration is a more highly pos- 
sible structure of the SDW on the basis of the obser- 
vation that there were no appreciable indications of a 
strain wave with 2Q and there was no third-harmonic 
component. Following that work, most of the ground- 
state calculations for 7-Fe were concentrated on finding 
a wave vector which minimizes the energy within the 
1Q helical structure. It should be emphasized, however, 
that the 3Q helical MSDW with Q = (0.1,0, l)2ir/a, 
(1, 0.1, 0)2n/a, and (0, l,0.1)27r/a is also consistent with 
the experimental results. This is because the neutron 
diffraction analysis cannot distinguish between the 1Q 
and 3Q states 40 when the crystal structure of the 7- 
Fe precipitates is properly cubic and the distribution 
of domains is isotropic. The MD approach presented by 
Kakehashi and co-workers 3 1-33 can predict the ground 
state without assuming the magnetic structure at the 
beginning. The resolution for a wave vector of magnetic 
structure in the MD calculations, however, is S > 0.2 
(in units of 2n/a) at the present stage; therefore it 
cannot reproduce the MSDW having the observed frac- 
tion (5 = 0.1). 5 One needs more accurate ground-state 
electronic-structure calculations to allow for the possibil- 
ity of the 3Q helical MSDW with the experimental wave 
vector. 



Regarding the consistency between theory and exper- 
iment, it should also be noted that although the experi- 
mentally suggested magnetic structure is a helical SDW, 5 
one should not exclude the possibility of the commensu- 
rate 28 ' 33 and linear 33 MSDWs which were found in the 
ground-state calculations for 7-Fe, because of the strong 
dependence of the magnetism of 7-Fe on the volume and 
strain. Experimentally, the SDW of cubic 7-Fe is found in 
a narrow range of lattice constants approximately equal 
to that of Cu, and the volume dependence of the mag- 
netic structure of cubic 7-Fe has not been investigated. 
It is also noted that the possibility of a small lattice dis- 
tortion is suggested at the onset of the IQ SDW in the 
7-Fe precipitates in Cu, 45 which might change the stable 
structure of 7-Fe. 

In the present phcnomcnological analysis, we focussed 
upon the magnetism of fee transition metals, specifically, 
that of cubic 7-Fe. Because the spin-orbit coupling effects 
are small in these systems, we neglected the anisotropic 
terms C(l,l',l",l"') in free energy (2.2). In order to ex- 
amine the effect of anisotropy, we also calculated the 
magnetic phase diagrams including the anisotropic terms 
in the free energy. We found two main effects. First, in- 
clusion of the anisotropic term partially removes the de- 
generacy of each SDW state; the most stable states be- 
come the SDWs having the longitudinal and transverse 
polarizations with respect to the x, y, and z axes, or the 
superposition of the three longitudinally (transversely) 
polarized states. Note that the longitudinal and trans- 
verse SDWs are still degenerate there relative to each 
other. Second, the phase boundary between SDW states 
is subject to a small displacement due to the anisotropy, 
but the global features of the magnetic phase diagrams 
in Figs. 1-4 remain qualitatively the same as long as the 
anisotropic terms are sufficiently small; the conclusions 
of the present work are not changed by considering the 
anisotropic terms. 
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Appendix A: Derivation of Free Energy (2.2) 

In this appendix, we derive the expression of real-space 
free energy (2.2) using the symmetry argument for the 
cubic system. We start from free energy expansion (2.1): 

N x,y,z 

f = ip Y Y aa p( l > l ') m i« m i>p 

1,1' a,0 

N x,y,z 

+ ]y4 Y Y b a /3ys{l,l',l",l'") 
1,1', I", I'" a,0,f,6 

x m to m i / /3 TO Z // 7 m ; /// 5 . (A-l) 

For the fee system, free energy (A-l) must be invari- 
ant with respect to ir/2 rotations of magnetic moments 
about the x, y, and z axes. It follows that the expansion 
coefficients satisfy the following relations: 



12 J. Phys. Soc. Jpn. Full Paper T. Uchida and Y. Kakchashi 



flaafi, 0' s are the same for a — x,y, z, (A-2) 



(y,z)(z,x)(x,y) 

aTYli' oJTliii fiTYli'ii ft 



o a/J (Z,O = fora^/3 (a, /3 = a;, y,z), (A-3) +m I/JWa m,». a )] . (All) 

b aaaa (l, I', I", l"'ys are the same for a = x,y, z, (A-4) Thug wg reach gq _ , 22 ^ 

&P( aa ^)(I, Z', Z", Z'") =0 for a 7^ /3 (a, /3 = x, y, z), In order to obtain the expression for an isotropic sys- 

(A-5) tern which is invariant under any rotation of magnetic 

i /f //t moments, it is sufBcient to require that free energy (2.2) 

baa/3 (Z, 1,1,1 ) s are the same ]-, e invariant with respect to the 7r/4 rotation of magnetic 

for (a,0) = (y,z)(z,y)(z,x)(x,z)(x,y)(y,x), (A-6) moments about the z axis: 
baf3f3a(l, I', I" , l'")'s are the same 

for [a,f3) = (y,z)(z,y)(z,x)(x,z)(x,y)(y,x), (A-7) f({m lx },{m ly },{m lz }) 

b a0 a {l, I', I", H's are the same = f({^=(m lx - m ly )}, {^(™ix + m ly )}, {m lz }). 

for (a,(3) = {y,z){z,y){z,x)(x,z){x,y){y,x), (A- 8) ' " ( A . 12 ) 

bp( aa fr)(l,l',l",l"') = This yie i ds the relation C{l,V ,1" ,V") = 

for a 7^ /3 7^ 7 (a, /3,-f = x,y, z). (A-9) b yyzz (l, I, I", I'") + b yzzy (l, I ,1 , 1) + b yzyz (l, I", I ,1 ) — 

b xxxx (l,l' ,1" ,1"') — 0. The resulting free energy satisfies 

Here, the subscripts P(aaa0) in eq. (A5) and P(aai3j) thc rotational invariance and thus describes the free 

in cq. (A-9) denote permutation of indices in thc parcn- energy for thfi isotropic systcm 
theses. By virtue of thc relations (A-2)-(A-9), free energy 

(A-l) is written as Appendix B: Expression of Free Energy with 

Three Wave Vectors 

Xi y jZ The phenomenological free energy for fee crystals is 

f = J_ ^ a xx (l, I') ]T m iami'a S iven b y ec l- ( 2 - 4 ) in § 2 for the general q vectors. In this 

u> a appendix we derive the expression of the free energy for 

SDWs consisting of three wave vectors, Q l , Q 2 , and Q 3 . 



- y 



iV 4 



u xxxx 

(I, I', I", I'") ^2 miami> a mi"ami>> 



WW" 

(y,z)(z,x)(x,y) 



+ b yyzz {l,l',l",l"') £ (miami'ami"i3mi"'i3 



(«,/?) 



B.l Commensurate case 

We consider here the SDW structures having commen- 
surate wave vectors 

Q l = (1,0, 0)27r/a, Q 2 = (0, 1, 0)27r/a, 



+ m w m ll pm l „ a m lllla ) Q 3 = (0, 0, l)27r/a. (B-1) 

+ b y zzy(l,l',l",l"') E (rruam^rm^m^a The magnetic moment vectors m(0i), ™(Q 2 ), and 



(«,/3) 

+ mipmi>ami"ami'"i3) 



m(Q 3 ) are real and orthogonal to each other: 



m(Q 2 )-m(Q 3 ) = m(Q 3 )-m(Q 1 ) = m(Q 1 )-m(Q 2 ) = 0. 

(y,z)(z,x)(x,y) (B-2) 

+ b yzyz (l,l',l",l"') £ (mi 

a m v /3 m i" a m i"' j3 The free energy having three wave vectors can be written 
with the help of eq. (2.4). Using relations (B-l) and (B-2), 
we reach 



+mipmv a m V ipm V u a )\ . (A-10) 

By grouping the same fourth-order terms of {to/ q } in the 
summations with respect to site indices, we have 



/ co = 4^[A(Q,)|m(Q l )| 5 



»=i 



1 — ^ +(i? 1 (Q 4 ) + B 2 (Q 4 ))|m(Q^I 4 l 



mialTll'a 



a (2,3)(3,1)(1,2) 

+ E B 3 (Q t ,Q 3 )\rnm\ 2 \rn(Q 3 )\ 2 . (B-3) 

u xxxx 

771 / Q 111 1 1 a 111 i" a 771 1'" a 

The coefficients B^Q,), B 2 (Qi), and B 3 (Qi, Qj) are de- 
+{Z>j/j/zz(Z, Z , I", l"')+b yzzy (l, I", I , l')+b yzyz (l, I", I , Z ")} fined as follows: 



iV 4 

IJ'V'l" 
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This is the free energy for the commensurate case given 
by eq. (3.10) in §3. 

B.2 Incommensurate case 

In the incommensurate case, we consider three equiva- 
lent wave vectors, Q 1 , Q 2l and Q 3 , which satisfy incom- 
mensurate conditions (4.41)-(4.43) in §4. Using expres- 
sion (2.4), the free energy having three incommensurate 
wave vectors can be written as follows: 



JMOi) = E B({Qi, -QJ. {Q it -QJ) (i = 1, 2, 3), 

(B-4) 

S 2 (Qj) = B(Qi, Qi, Qi, QA 

+ B(-Q i ,-Q i ,-Qi,-Q i ) 

+E[ B «(^)> 

p 

+ B({(Q i ,Q i ),{Q i ,-Q i }}) 
+B({(-Q.,-Q.),{Q.,-Q.}})] (i= 1,2,3), (B-5) 

B 3 (Q l ,Q J ) = E [^({{Qi.-Qi}.^)-^}}) 

p 

+ B({(Q j ,Q i ),{Q j ,-Q j }}) 
+ i?({(-Q.,-Q.),{Q J .,-Q J .}}) 

+ B({{Q i ,-Q i },(Q J .,Q J .)}) 
+ B({{Q.,-Q.},(-g.,-g.)}) 

+ B({(Q i ,Q i ),(Q j ,Q j )}) 
+ B({(-Q i ,-Q i ),(-Q j ,-Q j )}) 
+ B({(Q i ,Q i ),(-Q j ,-Q j )}) 
+B({(-Q i ,-Q i ),(Q j ,Q j )})~ 

((m) = (2,3)(3,1)(1,2)). (B-6) 

In eqs. (B-4)-(B-6) and below, X) p denotes summations #i(Qi, Qj) = E S (UQ 4 , -QJ> {Qj, ~Qj}}) 
with respect to all the permutations of elements in the p 
curly brackets. Wherever a parenthesis or curly bracket 
appears, it is regarded as one element when counting 
permutation. 

Since the wave vectors Q 1 , Q 2 , and Q 3 are all equiv- 
alent on the fee lattice, it follows that 



/ te = E^QJMQi)! 2 + Bi(Q i )|m(Q i )| 4 

i=l 

+ i? 2 (Q l )m 2 (Q t )m* 2 (Q 4 )] 

(2,3)(3,1)(1,2) 

+ E [Si(Q i ,Q j )|m(g i )| 2 |m(Q j )| 2 

+ B2(Q 4 ,Q,)|m(Q t )-m(Q J )| 2 

+ B 3 (Q i ,Q J -)MQi)-"»*(Q J -)| 2 ]- (B-12) 

Here, the coefficients B^QA, B 2 (Qj) (i = 1,2,3), 
B 1 (Q i ,Q j ), B 2 (Q i ,Q j ), and B 3 (Q t ,Q 3 ) = 
(2, 3) (3, 1)(1, 2)) are defined as follows: 

B 1 {Q i ) = E B({Q 4 , -QJ, {Qi, QJ) (i - 1, 2, 3), 

(B-13) 

B 2 (Q l ) ee J2 B({(Qi, Qi), ("Qi, -Qi)}) (< = 1, 2, 3), 

(B-14) 



((m) = (2,3)(3,1)(1,2)), 

B 2{Qi, Qj) = E Qj}> {-Qi. -Qi») 



i Q = 44(QJ = 4A(Q 2 ) = 4A(Q 3 ), 
S 1Q = Si(Qi) = Si(Q 2 ) = Bi(0s). 
B 2Q = B 2 (Qi) - B 2 (Q 2 ) - B 2 (0 3 ). 



(B-7) 
(B-8) 
(B-9) 



((M) = (2,3)(3,l)(l,2)), 
B 3 (Qi, Qj) = E B («Q" -Qj}> {"Qf Qj}» 



(B-15) 



(B-16) 



Biqq = B 3 (Q 2 , Q 3 ) = B 3 (Q 3 , Qi) = B 3 (Qi, Q 2 )- 

(B-10) 

With the use of relations (B-7)-(B-10), free energy (B-3) 
is expressed as 



/co = E^QMQ,)! 2 + ( B iQ + BiQ)\™{Qi)\ 4 



((i,j) = (2,3)(3,l)(l,2)). (B-17) 

Since the wave vectors Q 1; Q 2 , and Q 3 are all equivalent 
in the fee lattice, it follows that 



A Q ee 2A(Q 1 ) = 2A(Q 2 ) = 2A(Q 3 ), 
B 1Q ee BiiQJ = Bi(Q 2 ) = Bi(Q 3 ), 
B 2Q ee B 2 (Q 1 ) = 5 2 (Q 2 ) = B 2 (Q 3 ), 



(B-18) 
(B-19) 
(B-20) 



+ 



(2,3)(3,1)(1,2) 

E 



Bi QQ |m(Q 4 )| 2 |m(Q,)| 2 - (B-ll) 



Bi QQ ee Bt(Q 2 , Q 3 ) = B^Qs, QJ = B 1 (Q 1 ,Q 2 ), 

(B-21) 

S 2QQ = B 2 (Q 2 , Q 3 ) - 5 2 (Q 3 , = B 2 (Qi, Q 2 ), 

(B-22) 
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B 3QQ = B 3 (Q 2 ,Q 3 ) = B 3 (Q 3 , = B 3 (Q U Q 2 ). 

(B-23) 

With the use of relations (B-18)-(B-23), free energy 
(B-12) is expressed as 



/ ic = ^[A Q |m(Q l )| 2 + B 1Q |m(Q J )| 4 
+ B 2Q m 2 (Q i )m* 2 (Q i )} 

(2,3)(3,1)(1,2) 

+ X! [SiqqI^QJIV^OI 2 

+ B 2QQ \m(Q l )-m(Q J )\ 2 

+ B 3QQ \m{Q l )-m*{Q J )\ 2 ]. (B-24) 

This is the free energy for the incommensurate case, 
eq. (4.44), in §4. 
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